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We use the asymmetric simple exclusion process for describing vehicular traffic flow at the in- 
tersection of two streets. No traffic lights control the traffic flow. The approaching cars to the 
intersection point yield to each other to avoid collision. This yielding dynamics is model by imple- 
menting exclusion process to the intersection point of the two streets. Closed boundary condition 
is applied to the streets. We utilize both mean-field approach and extensive simulations to find the 
model characteristics. In particular, we obtain the fundamental diagrams and show that the effect 
of interaction between chains can be regarded as a dynamic impurity at the intersection point. 
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I. INTRODUCTION 



Modelling a vast variety of non equilibrium phenom- 
ena has constituted the subject of intensive research 
by statistical physicists [2, H[. ^ n particular, vehicu- 
lar dynamics has been one of these fascinating issues 
[3> 0) 111- While the existing results in highway traffic 
needs further manipulations in order to find direct appli- 
cations, researches on city traffic [H, 0, HI seem to have 
more feasibility in practical applications. Recently, no- 
table attention have paid to controlling traffic flow at in- 
tersections and other d esig nations such as roundabouts 
d 01, EE2, E 0, B H||. In this respect, we intend 
to study another aspect of traffic flow at intersections. 
In principle, the vehicular flow at an intersection can be 
controlled via two schemes. In the first scheme the traffic 
is controlled without traffic lights. In the second scheme, 
signalized traffic lights control the flow. In the former 
scheme, approaching car to the intersection yields to the 
traffic in its perpendicular direction by adjusting its ve- 
locity to a safe value to avoid collision. The basic ques- 
tion is that under what circumstances the intersection 
should be controlled by traffic lights? In order to cap- 
ture the basic features of this problem, we construct a 
simple stochastic model. The vehicular dynamics is rep- 
resented by asymmetric simple exclusion process (ASEP) 
[HI, G3, EH- The intersection point is the place where 
two chains representing the streets interact with each 
other. It is a well-established fact that a single static 
impurity can strongly affect the characteristics of ASEP 
both in closed 0, [2(| and open boundary condition [2l[ . 
In addition, the characteristics of ASEP in the presence 
of moving impurities has been studied and shown to ex- 
hibit disorder- induced phase transitions [13] • Besides rel- 
evance to traffic flow, the investigation of ASEP in the 
presence of small amount disorder has recently revealed 
the existence of novel aspects of the interplay of disor- 
der and drive [H [H HI, HI] . In our model, the effect 
of the perpendicular chain can be interpreted as a single 
dynamic site-wise disorder which to our knowledge has 
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FIG. 1: Two perpendicular closed chains representing the in- 
tersection of two uni directional traffic flows. 



not been investigated. 



II. 



DESCRIPTION OF THE PROBLEM 



Consider two perpendicular one dimensional closed 
chains each having L sites (L is even) . The chains repre- 
sent urban roads accommodating unidirectional vehicular 
traffic flow. They cross each other at the crossing sites 
i\ = «2 = 2" on the first and the second chain respec- 
tively. With no loss of generality we take the direction of 
traffic flow in the first chain from south to north and in 
the second chain from east to west (see Fig.l for illustra- 
tion). Each site of the chain is either vacant or can hold 
at most one car (hereafter interchangeably called par- 
ticle). We assign an integer valued occupation number 
rii (rrii) to site i of the first (second) chain respectively. 
In case the site is occupied by a particle, its occupation 
number is one and zero otherwise. 

The system configuration at each time t is character- 
ized by specifying the occupation numbers rii and rm 
(i = 1, • • ■ , L) in both chains. The system dynamics is 
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asymmetric simple exclusion process (ASEP). Although 
this choice is far from the realistic vehicular dynamics, 
we expect the general features of the problem namely the 
phase structure does not qualitative change. In addition, 
employing ASEP dynamics allows us to treat the problem 
more easily on analytical grounds. In reality, each driver 
yields to perpendicular traffic flow by appropriately ad- 
justing of her velocity. In our formulation, this caution- 
ary behaviour can be described by simple exclusion pro- 
cess. This process is defined as follows. During an in- 
finitesimal time dt each particle can stochastically hop to 
its forward neighbouring site provided the target site is 
empty. If the target site is already occupied by another 
particle, the attempted movement is rejected. At each 
time t, we define the particle approaching to the inter- 
section as the one occupying the previous site % = ^ — 1. 
We note that this particle's attempt is successful if both 
sites i\ = t£ and i 2 = 4 which is actually the common 
intersection site on both chains is empty. In terms of 
ASEP terminology, the sites 4 — 1 can be regarded as 
dynamic impurity sites for each chain. The hoping rate 
at all sites are scaled to unity except at i\ = 4 — 1 and 
%2 = -j — 1 which are 1 — rriL and I — til respectively. 



A. Mean field approach 



(l-rii)(l-mi). The latter expression is the probabil- 
ity that the middle sites in each chain are simultaneously 
empty. In the steady state, the left hand sides of mean- 
field equations become zero and we arrive at a set of 2L 
nonlinear equations. Even by employing the assumption 
of mean field, we are not able to solve these nonlinear 
algebraic equations. Therefore, we should resort to nu- 
merical methods. We now outline a numerical approach 
for solving the set of nonlinear equations. 

B. Numerical approach to mean field equations 

Our approach for solving the MF equations is based on 
the constant density scheme which has originally been in- 
troduced by Barma and Tripathy [UHU . In this scheme, 
we first fix the global densities in two chains at given val- 
ues pi and p2- Second, we assign initial density profiles 
ni[0], n 2 [0], ■ ■ ■ ,ni[0] and mi[0], to 2 [0], • • • , tol[0] to the 
first and second chain respectively. The constancy of 
global densities implies the following constraints on the 
initial profiles: 

n 1 [0]+n 2 [0] + --- + n L [0]=Lp 1 (4) 



Let us denote the mean density at site i and time t of 
the first and second chains by (m) and (to,) respectively. 
The master equation governing their time evolution can 
simply be written as follows: 



mi [0] + m 2 [0] + ■ • • + m L [0] = Lp 2 (5) 

We next evolve the site densities according the follow- 
ing discrete time updating rules: 



— (rij) = (nj_i(l - rii)) - (n;(l - n i+1 )) (1) 

In the above equation, the site index i covers all the 
lattice sites except i = 4 — 1 and « = The rate equa- 
tion for the second chain can be obtained by replacing m 
by rtii . For sites i = 4 — 1 and i = 4 , the rate equations 
have the following forms: 



^("f-i) = (ni^fi-^-i)) - (n£_ x (l-ri£ - mt)) 

(2) 



n l [t + l] = n l [t]+n^ l [t]{l-n l [t])-n l [t]{l-n l+l [t]) (6) 

Similar equations hold for m^. Note that in the above 
equations, i covers the whole chain except the sites 
2 = 4 — 1 and i = 4^. The above discrete time evolu- 
tion rules stem in time discretisation of the mean-field 
equations within Euler algorithm. In the special sites 
i = j — 1, j the interaction between two chains modifies 
the rate equations and we should take into account equa- 
tions (2,3) which gives rise to the following dynamical 
rules: 

+ 1] = rii-iM + riL_ 2 [t}(l - nL^t])- 



dt 



(til) = (riL^l-riL - mi))-(nt(l-tii +1 )) (3) 



(7) 



.[t + l}=n,[t}+n,_ i m-ndt})(l 



Similarly, the rate equations for tul _ x and ttil can be 
obtained by replacing m <-> n respectively. In order to 
proceed analytically we take into account the mean-field 

approximation. In this approximation, we replace the ul [t](l — nt , -Ji]) (8) 

two-point functions by the product of one-point func- 
tions and furthermore we replace the probability that The equations for m,L_ x and tul are simply obtained 
the middle site i = % is empty, i.e. (1 — rii - mt), by from the above equations via replacing m by n and vice 
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FIG. 2: First chain current Ji versus its density pi for various 
densities p2 in the second chain. L = 300. 



versa. Notice that our dynamical equations preserve the 
constancy of global densities. More concisely, we have 

m [t] + 712 [t] H 1- til [t] = Lpi and mi [t] + m% [t] H h 

rai[i] = Lp2 at each t. 

With an appropriate choice of initial condition, after 
iterating the above equations for many time steps the sys- 
tem is expected to reach a fixed point denoted by {n*} 
and {m*} in which further iteration does not change the 
densities. This solution can be considered as the solution 
of the mean-field equations. Note that in a acceptable 
solution, all the densities n\ , ■ ■ ■ ,n* L and m\ , ■ ■ • , m* L 
should lie between zero and one. The chains currents are 
thus obtained according to the relations J\ = n*(l—n* +1 ) 
and J 2 = m*(l — m* +1 ) in which i can be any site of 
the chains. It should be noted that the choice of initial 
densities is a crucial step. Only certain initial condition 
converges to the desired solution. In general, the long- 
time behaviour of the density profile turns out to be an 
oscillatory pattern. 



C. Monte Carlo simulation 

For obtaining a better insight, we have also executed 
extensive Monte Carlo simulations which are presented 
in this section. The chains sizes are equally taken as 
L\ = L2 = 300 and we averaged over 100 independent 
runs each of which with 10 5 time steps per site. After 
transients, two chains maintain steady-state currents de- 
noted by Ji and J 2 which are function of the global den- 
sities pi and p2- We kept the global density at a fixed 
value pi in the second chain and varied p\. Figure (2) 
exhibits the fundamental diagram of the first chain. 

The generic behaviour is reminiscent to ASEP with a 
single defective site [lj|[2(J[28j]. Intersection of two chains 
makes the intersection point appear as a dynamical de- 
fect. It is a well-known fact that a local defect can affect 
the system on a global scale 0, [2l[ . This has been con- 
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FIG. 3: Second chain current J2 versus density p\ of the first 
chain for various densities p2. 



firmed not only for simple exclusion process but also for 
cellular automata models, such as Nagel-Schreckenberg 
[29l | , describing vehicular traffic flow [j, [3(1 [HJ . Analo- 
gous to static defects, in our case of dynamical impurity 
we observe that the effect of the dynamic defect is to 
form a plateau region p 6 [p_ , p+] in which p± = 0.5 ± A 
and 2A is the extension of the plateau region in the fun- 
damental diagram. In the plateau region the current is 
independent of the global density. The larger the den- 
sity in the perpendicular chain is, the dynamic defect has 
larger strength. For higher p 2 , the plateau region is wider 
and correspondingly the current value is more reduced. 
The notable point is that the flow capacity in the first 
chain persist to large decrease up to considerably large 
density p 2 ~ 0.5 in the second chain. This marks the fact 
that conflicting flows of particles can, to a large extent, 
weakly affect each other. To shed some more light on this 
aspect, let us now consider the flow characteristics in the 
second chain. In figure (3) we sketch the behaviour of 
perpendicular fundamental diagram that is J2 versus p\ . 

As depicted, J 2 shows a phase transition at the criti- 
cal density p±. c = p+. Before pi. c , the current exhibits 
smoothly decreasing behaviour. The nature of decrease 
depends on the value of p2- For p2 < 0.3 or p2 > 0.7 
J2 is almost constant and is obtained from the single 
chain relation J 2 = P2O- ~ Pi)- However, in the inter- 
val 0.3 < P2 < 0.7 J 2 shows a complex behaviour as 
is shown in figure (3). It first increases up to a small P2 
then smoothly diminishes until it reaches to a plateau re- 
gion. The reason is due to interaction between two chains 
which induces correlations between them. This modifies 
the value of J2 from the single chain value P2(l — p2)- The 
appearance of a maximum in J2 marks the point that a 
small density of cars in the first chain can even regulate 
the traffic, and hence enhance the flow, in the second 
chain. When the global density in the first chain exceeds 
the critical value, the perpendicular current exhibits a 
quasi linear decline. This corresponds to capacity break 
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FIG. 4: Density profile by MC simulation in the first chain: 
pi = 0.2 and pi — 0.4 (top) and pi = 0.6 and p2 = 0.8 
(bottom). 



down in the second chain. In terminology of vehicular 
traffic, if the density of the perpendicular chain goes be- 
yond a critical value, one should be warned that control- 
ling of the traffic via self-organised mechanism starts to 
fail and traffic lights signalisation is thereby prescribed. 



D. density profiles 

In order to improve our understanding, it would be 
useful to look at the behaviour of density profiles in both 
chains. Let us look at some typical density profiles before 
attempting to give a general remark. The following fig- 
ures, obtained by simulation, display the density profile 
in the first chain for two densities p\ — 0.2 and p\ = 0.6. 

In figure (4) top, the second chain affects the first 
chain's density profile only on a local scale. The reason 
is that the p\ is below the limit to be globally affected. 
On the other hand, if one increases p\ to a higher value 
e.g. 0.6 (figure (4) bottom), it turns out that density pro- 
file of the first chain gives rise to segregation between a 
high and a low density regions. The formation of phase- 
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FIG. 5: Density profile in the first chain obtained via MC 
(dashed) and MF (dotted), pi = 0.4 and p 2 = 0.8. 



segregated regime depends on the mutual values of global 
densities in both chains. The above observations have 
resemblance to the BML model of city traffic @ in the 
sense that below a given density, the conflicting flows do 
not affect each other much and the cars are not notably 
blocked by the other lane flow. Nevertheless, it should 
be mentioned that updating rules and road structures 
in the BML model are entirely different to our model's. 
Furthermore, the blocking mechanism in the BML model 
is not only due to exclusion principle but is also related 
to the cooperative motion of vehicles. In our model, the 
it is only the exclusion at the crossing point which gives 
rise to blocking. To support our density profiles find- 
ings, which have been obtained by Monte Carlo simula- 
tions, we have numerically solved the mean-field equation 
in the constant-density scheme described in section II. B. 
The results are satisfactory and in qualitative agreement 
with Monte Carlo simulations. As an example, in figure 

(5) we have sketched and compared the profile of density 
in chain one for given global densities in both chain both 
with Monte Carlo and numerical mean field approach. 

We observe that numerical mean field has qualitatively 
reproduced the phase-segregated behaviour. The differ- 
ence between high and low densities given by MC and 
MF are in fairly well agreement. However, there is a no- 
table difference in the vicinity of the domain wall. The 
prediction of MF is much sharper than that of MC. The 
more smooth transition from low to high density in MC is 
related to involve fluctuations which are not captured by 
the MF approach especially in low dimensions. In figure 

(6) we show density profiles in uniform phases. 

The predictions of MF and MC are reasonably close to 
each other. Some comments on the constant-density ap- 
proach seems unavoidable. In each iterative scheme, the 
stability of long time behaviour should be investigated. 
In our case, the choice of initial values of density profiles 
play an important role. We took their values within a 
small interval of size r\ around the corresponding global 
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FIG. 6: Density profile in the first chain obtained via MC and 
MF. pi = 0.2 and p 2 = 0.6. Number of iterations is 2000. 
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FIG. 8: Phase diagram of the model for the first chain, o 
shows the corresponding diagram for the second chain. 
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FIG. 7: Density profile in the first chain obtained via numer- 
ical MF. pi = 0.2 and p 2 = 0.6. 



densities. In figure (6) we set ?y = 0.01. We note that the 
long time behaviour of the system of equations give rise to 
an oscillatory pattern of profile. In numerical terms, the 
iteration becomes unstable for large times. However, in 
intermediate times, iterative method gives a reasonable 
answer. In figure (7) we have depicted the behaviour of 
the density profiles for a larger iteration. We have varied 
the the number of iterations M for 2000, 5000 and 10000 
(from bottom to top) correspondingly. For clarity, den- 
sities are shifted upwards for M = 5000 and M = 10000. 

We now turn into the general behaviour of density pro- 
files. It would be a natural question to ask under what 
circumstance yielding leads to traffic jam formation in 
the first chain before the crossing point. To this end, by 
Monte Carlo simulations we have systematically surveyed 
the entire phase space p\ — pi in the grids of 5p = 0.05 
. Two phases are identified: uniform (homogeneous) and 
phase-segregated. In the uniform phase, the interaction 
between two chains has only a local effect on profile of 



the first chain whereas in the segregated phase, a domain 
wall separates a low density and a high density region. 
Figure (8) exhibits the phase structure of the problem. 
By symmetry, one obtains the phase structure of the sec- 
ond chain via replacing 1 «-> 2. 



III. MEAN-FIELD PHASE DIAGRAM: 
ANALYTICAL APPROACH 

Here we try to outline an analytical approach to obtain 
the model phase diagram. This approach is based on the 
concept of density phase segregation proposed in [n| [2(| • 
Suppose for given densities p± and P2 one has density 
phase segregation in both chains. We denote the high 
and low densities in the first chain by by ph and pi and 
in the second chain by ^ and respectively. Also let us 
denote the relative length of low and high density regions 
in the first chain by a; and ah and in the second chain 
by bi and bh- We denote the probability of occupation of 
the crossing site by a car from the first (second) chain by 
T\ (t2)- Therefore the probability that the crossing site 
will be empty is 1 — r± — r 2 ■ We can write the following 
mean field equations: 



Ji = Pi(l - Pi); Ja = 6(1 - CO- 



(9) 



Jx=p h {l-Ph)\ Ja = &(!-&) (io) 



J 1 = p h (l-n-r 2 y, J 2 =a(l-n-r 2 ) (11) 



Ji=n(i-pO; ^2 = r 2 (i-6) (12) 
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FIG. 9: Phase diagram of the model. 



a t +a h = 1; 



(13) 



piai + phdh = pv, &bi + £hbh = P2 (14) 
The above equations are easily solved and we find: 



In the interior region, one has density phase segregation 
in both streets. We note that middle square region is 
in well qualitative agreement with the prediction of MC 
simulations ( intersection of two curves in fig8 ) . 



IV. 



SUMMARY AND CONCLUDING 
REMARKS 



We have investigated the characteristics of two con- 
flicting traffic flows within the framework of asymmetric 
simple exclusion process. Two perpendicular chains in- 
teract with each other via the intersection point. Using 
Monte Carlo simulations and numerics, we have obtained 
the dependence of each chain current on its own and on 
its perpendicular chain global density. It is verified that 
the chains can maintain large currents up to rather a 
high density. Interaction of two chains can effectively be 
considered as a dynamic impurity. For some values of 
global densities in the chains, the interaction of chains 
leads to formation of high density region behind the in- 
tersection point which is segregated from a low density 
region afterwards. By a systematic scanning of the phase 
space, we have obtained the structure of model phase di- 
agram. Two phases of jamming (density segregated) and 
regulated (uniform density) flows are identified. 



1 



n = r 2 



pi = 6 



l 



Ph = 6i 



Jl = J 2 



15) 



One also finds: ai = 2 — 3pi and bi = 2 — 3p 2 - The 
conditions < a; < 1 and < bi < 1 imply: -| < pi, P2 < 
|. In figure (9), we have plotted the phase diagram. 
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